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I - INTRODUCTION 



Since the pionnering work of Anderson [1], the locahzation problem in one di- 
mension has been investigated extensively and is by now pretty well understood [2]. 
However, despite the large number of rigorous results, there are very few solvable mod- 
els in the continuum for which one can compute exactly the density of states and the 
localization length. The purpose of this work is to present some new results for such a 
model which is defined by the one dimensional Schrodinger Hamiltonian 

H=-^ + cj>\x) + cj>'{x) (1.1) 

This Hamiltonian, which was introduced by E. Witten as a toy model of supersymmetric 
quantum mechanics [3], has stimulated a number of interesting developments in the 
context of quantvmi mechanics. In particular it has provided exact solutions of the 
Schrodinger equation for a class of so called shape invariant potential [4]. It has also 
inspired a new method of semi-classical quantization [5,6]. In the context of disordered 
systems, i.e. when the superpotential (f){x) is considered as random, one of the most 
interesting features of this model comes from its relation to the problem of classical 
diffusion of a particle in a one dimensional random medium. This correspondence [7], 
which has been exploited in great detail in [8], is based on the observation that the 
imaginary time Schrodinger equation 

^ + Hi; = (1.2) 

can be cast into the Fokker-Planck equation 



dt dx \dx 
through the transformation 

Alternatively, the diffusion process is described by the Langevin equation 



x{t) = 2(j){x) + r]{t) 
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(1.4) 



where r]{t) is a gaussian white noise such that 



(1.5) 



This correspondence enables one to express the transition probabihty P{x t \y 0) in 
terms of the eigenstates of H. In the disordered case, namely when 0(a;) is random, the 
probability of returning to the starting site, averaged over the disorder, becomes 



^0 

The long time behaviour of P is thus related to the low energy behaviour of the density 
of states p{E). 

In the case of a white noise 

< (j){x) >— jia 

(1.7) 

< (j){x)(j){x') > - < (j)'^{x) >= aS{x - x') 

both the localization length and the density of states have been computed exactly [8,9] 

and display unusual behaviours at low energy. Denoting by N(E) the integrated density 

of states, one obtains N(E) oc — ^ — for u = , and N(E) oc for u > 0. 

E^oln E -E^o ^ 

This implies that < P{xt | x 0) > will present an algebraic tail at large time which 
reflects the anomalous behaviour of the position x{t). Since the anomalous behaviour 
of the diflFusion can be predicted qualitatively [10], the above correspondence can be 
used to understand physically the origin of the anomalous behaviour of the localization 
problem (see also [11,12]). 

However a direct interpretation of these results without appealing to such a cor- 
respondence has until now been missing. One of the purposes of this paper is to 
fill this gap. In our analysis, a crucial role will be played by the zero energy states 
ipo{x) = exp (j)(t)dt and ipi{x) = ipo{x) .^f x which are exactly known for any re- 
alization of the disorder. This very unusual feature in the context of disordered systems 
can of course be traced back to the supersymmetric structure of the Hamiltonian. 

In order to display this structure it is convenient to introduce the pair of Hamil- 
tonians 




(1.6) 



H± = 



dx"^ 



+ 02(a;)±0' 



(1.8) 
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They can be rewritten in the factorized form 

H+ = Q^Q and H_ = QQ^ , where Q = --^ + 6(x) (1.9) 

ax 

This imphes and H_ have the same spectrum for £■ > [3]. In the presence of 
disorder they are characterized by the same Lyapunov exponent and density of states. 
We will see that the possibility to treat them on the same footing will significantly 
simplify our analysis. 

In all this work, the disorder is modeled by assuming that the potential {^(a;)} 
is described by an ensemble of rectangular barriers with alternating heights 0o and 0i 
of random length. Most of our work deals with the case where {(f>{x)} is a random 
telegraph process [13]. The lengths of the intervals over which (j){x) is constant are 
therefore distributed according to an exponential law pi{£) = 9{£) niexp{—ni£) for 
i = 0, 1. Besides its mathematical interest, this form of disorder is motivated by the 
fact that it can model certain quasi-one-dimensional structures with a piecewise constant 
order parameter [14]. Another advantage is that this type of disorder is less singular 
than a pure white noise process for which arbitrary jumps of the potential can occur. 
An interesting extension of this model corresponds to the case where the lengths of 
the barriers are chosen according to a broad distribution which behaves for large £ as 

£-(l+a)_ 

This paper is organized as follows. In part II we present a method to compute the 
density of states which is an adaptation of the one discussed by Benderskii and Pastur 
[15]. We then apply this method to the case where the disorder is a random telegraph 
process. Exact analytical expressions of the density of states and localization length are 
derived. This is followed by a careful analysis of the limiting behaviour for £^ ^ 0. A 
physical interpretation of these results is given in part IV. Numerical simulations are 
presented in part V. Then the case where p{£) is a broad distribution is considered. 
Remarkably enough, both the density of states and localization length can be computed 
exactly when a < 1 and = — 0o- Finally, in part VI we apply some of these results 
to the study of anomalous diffusion and discuss various applications. 
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II- THE PHASE FORMALISM FOR THE SUPERSYMMETRIC 
HAMILTONIANS H± 



The phase formahsm is based on the foUowing property of one-dimensional Hamil- 
tonians : the integrated density of states per unit length N[E) and the Lyapunov 
exponent ^{E) are directly related to properties of the solution '^e of the corresponding 
Schrodinger equation H'^e = E'^e with a given logarithmic derivative at one point [16]. 
The number N{E) of states of energy lower than E, per unit length, is equal to the 
number of nodes of the wave function i/'e per unit length, and the Lyapunov exponent 
7(£') measures the exponential growth rate of the envelop of the wave function i^^e- 

For any spatially homogenous disordered potential with short range correlations, 
N{E) and ^{E) are self-averaging quantities. This means that they are the same for all 
realizations of the disorder with probability one. 

A convenient way to implement the phase formalism for the Schrodinger equations 



The phase 9e{x) contains all the informations about the oscillations of the wave func- 
tions (ue, ve) and the modulus pe{x) represents the envelope of these wave functions. 
The new dynamical variables 9e and pE will determine the value of N{E) and ^{E) 
respectively. 

N{E) is equal to the number of nodes oi ue (or ve) per unit length and can thus 
be written directly in terms of 9e{x) 



Note that [9e{L) — ^s(O)] is the total phase accumulated on the interval [0, L] and is 
therefore not an angle defined modulo 27r. ^{E) is the exponential growth rate of the 
envelope pE and reads therefore 




(2.1) 



is to introduce "polar" variables 



ue{x) = pe{x) cos9e{x) 
ve{x) = Pe{x) sin0£;(x) 




(2.2) 




(2.3) 
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Now in order to write the dynamical equations for 6e and pe we have to separate the 
cases E > and < 0. From now on we will set k = C\E\. 



1 — case E = +k^ 

In this case, the factorized forms (1.9) allow to rewrite the system (2.1) as two 
coupled first order equations Qu = kv and Q'^v = ku, which explicitly read 



Note that the equation for Be does not involve pe- 

Let us now present the direct resolution of these equations in {9e, Pe) for the 
special cases E = and £^ — > oo, before we turn to the general method to calculate 

N{E) and -f{E) for aU E > 0. 

a) case E = 

For k = 0, the dynamical equations (2.5) read 



The phase 9{x) cannot grow more than 7r/2 on the interval [0, L] since, for any {(p{x)}, 
the velocity vanishes whenever ^ = [modulo tt/2\. Therefore (2.2) yields immedi- 
ately A'"(0) = in accordance with the positivity properties of H±. 

To get 7(0), we can integrate the equation for d{x) 




(2.4) 



This gives the following dynamical equations for the "polar" variables 




(2.5) 




(2.6) 



tan 6'(L) -2 ("^ 4>{x)dx 

= e Jo 

tan^(O) 



(2.7) 
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and rewrite the equation for p as 

d\n.p dO cos 26{x) 



dx dx sin 29{x) 

The integration gives 

sin 2^(0) 



(2.8) 



7(0) - lim_ ( \ ) = ^l^rn^^ m 



L— >oo 



sin2^(L) 



(2.9) 



We now use (2.7); in the hmit L — > oo, the expression (f){x)dx is simply the mean 
value of {0(x)}, that we will note Fq =< > from now on. We obtain finally that the 
Lyapunov exponent at zero energy is equal to the absolute value of the average Fq of 
{cPix)} 

7(0) = |Fo| (2.10) 

This simple result can also be recovered by starting from the two linearly inde- 
pendent solutions of H+ip = 

, / N + r <t){t)dt 

ipo[x) = e Jo ' 



The asymptotic behaviour for large x 

/ (l){t)dt ~ Fox[l + o{l)] 

Jo a;^oo 



(2.11) 



allows to recover (2.10). 



b) Limit E ^ oo 



dO 

At high energy, (2.5) can be approximated by — ~ A; = CE and N{E) behaves 

dx fc— >oo 

CE 

therefore asymptotically as the free integrated density of states N{E) 



E—>oo TT 



At this order of approximation, the Lyapunov exponent vanishes 
^{E) = lim < ^ / dx(l){x) cos[26'(x)] > 

L^oo L Jq 

1 

~ lim < (f) > — dxcos{2kx) = 

— > oo L— >00 L Jq 



oo L— >oo 

7 



To find the asymptotic behaviour of ^{E) we must calculate the next order for 0{x) 

e{x) ^ kx- I (j){x') sm{2kx')dx' (2.12) 

^=^0° Jo 

This allows to relate the high energy behaviour of ■jiE) to the two point correlation 
function of the process {(j){x)} [16] 

1 r°° 

-f{E) ~ - / dx cos(2/ca;) [ < (p{x)(p{0) > - < >^ 1 (2.13) 

E^oo 2 Jo 

This Fourier transform relation shows that the high energy decay of the Lyapunov 
exponent is directly related to the regularity of the random process {(f){x)}. 

For the white noise case < (j){x)(j){x') >=< (f) >^ +a5{x — x') we recover the result 

[8] 

-^WN^E) f (2.14) 

This behaviour at high energy is rather pathological and very particular to white noise. 
Indeed for an exponentially correlated noise 

< cj){x)(j){x') >=< </. >2 +^e-^l^-^'l (2.15) 

which can be viewed as a regularization of the white noise case as long as A is finite, we 
obtain a vanishing Lyapunov exponent 



c) General method to calculate N{E) and ^{E) for all £^ > 



We have seen that the dynamical equations (2.5) could be integrated directly for 
the particular values E = ^ and — > oo in order to get NiE) and 7(-E). This is 
obviously not possible in general. We therefore need a more powerful method [15]. We 

— to rewrite the number Ml{E) of nodes 

Xi zero of f '"^ ^^^'^ 

of cosdE{x) in the interval [0, L] 

MhiE) = i dx 5{x-Xi) = k I dx 5[cosd{x)] 

Jo 1 ,■ Jo 

XI solution 

of cos e(a:j)=0 

j-L 

= k dxJ2 ^[^(^) - 2 ~ <^2.17) 



The integrated density of states N{E) then reads [15] 



N{E) = lim 



( 



L 



) 




(2.18) 



where Peq{d) is the stationary distribution of the reduced phase 9 defined modulo tt. 
Similarly the Lyapunov exponent can be rewritten as [17] 



where the brackets < ■ ■ ■ > denote the mean value taken over the joint probability of 
the process (f) and of the reduced phase 9. We will use these results in section III to 
calculate explicitly N{E) and ^{E) for a specific model of disorder. 

2 — case E = 

We have already seen that A'"(0) = so that N{E) vanishes identically for E e 
] — oo, 0]. It is nevertheless interesting to study the Lyapunov exponent 7(-E') in this 
unphysical region since it is relevant for the problem of classical diffusion in the random 
force field {0(a;)}. 

For E < 0, the factorized forms (1.9) allow to rewrite the system (2.1) as two 
coupled first order equations Qu = kv and Q^v = —ku, which explicitely read 




(2.19) 




(2.20) 



where k — C—E. 



The dynamical equations for the polar coordinates then read 




(2.21) 



Note the differences with (2.4) and (2.5) respectively. 
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We recover immediately that N{E) = for < by noting that the accumulated 
phase ${x) will be trapped in an interval of length < ^. Indeed for any {(/){x)} the 
velocity is positive whenever 9 — ^ [modulo tt], and negative whenever — 
[modulo tt] . 

Still we need to study the stationary distribution Peq{0) in order to calculate ^{E) 
along the lines explained in II 1 c). 

Ill - APPLICATION TO A MODEL OF RECTANGULAR BARRIERS 
OF RANDOM LENGTHS 



1 — Description of the model 



Let us now consider a model [15] where (t){x) takes alternatively two values 0o 
and 01 on intervals whose lengths are positive independent random variables (Fig.l), 
distributed according to the following probability densities respectively 

/o(/)=^(Orioe-"^ ^ ^ 

3.1 

h{h) = e{h) m e-^^^ 

This choice for /o and /i is in fact the only one that makes the process {0(x)} Markovian. 
This property enables us to write differential equations for the probability Pq{x) to have 
0(a;) = (^0 and the probability p\{x) = 1 — Pq{x) to have 0(a;) = 0i. 

dpQ 




-noPo + nipi = ni - (no + ni)po 
-nipi + noPo = no - (no + ni)pi 

The corresponding stationary solutions are simply 

ni < / > 



(3.2) 




<l> + < h> 

< h> 
<l> + < h> 



(3.3) 



where < / >= dl lfo{l) = — is the mean length of intervals {(j){x) = (po} 

^ no 

and < h >= dh hfi{h) = — is the mean length of intervals {(f){x) = 0i}. 

ni 
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The mean value Fq and the two point correlation function G{x) of the process 
{(j){x)} read 



no + ni no + rii 

Gix) <P{x)m > - < >2= ^"^^ (01 - 0o)'e-("°+"^)l-l 

(no + ni)^ 

The case 0i = — 0o where 

^(a;) = (0§ - F2)e-("o+ni)|x| 
tends to the white noise process 

< >= Fo 



(3.4) 



(3.5) 



in the limit 



00 - 


-> oo 




no - 


->■ oo 


where 


ni - 


-> oo 





< 0(x)0(O) >= F(f + c75(x) 



and remain constants. (3-6) 



(7 = 2 



It will be also useful to define /i = ^ = 



no+n-L 

24>o 



We will now see that this model is exactly soluble, even though it contains corre- 
lations, because it is a Markovian process. 

Note that this two-step process can be easily generalized into a multi-step process 

[18]. 

2 — Application of the Phase Formalism for E > 

As we have already stressed, {0(a;)} is a Markovian process. Since the accumulated 
phase 6{x) evolves for F > according to ^ — k — (j){x) sin26'(a;), we see that the pair 
{0(a;), 0{x)} forms a two-dimensional Markov process. 

Let us define 

e{x) e[e,e + de] 

Pq{9,x) dO = the probability to have \ and 

0(x) = +00 

e{x) e [6*, 6* + de] 

Pi{9,x) dO = the probability to have and 



0(x) = +01 
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We can write two coupled master equations for Pq and Pi 

- 00 sin 2e)Po] - noPo + mh 

(3.7) 







dx 


de 


dPi 


d 




de 



We now express the probability distributions Pq{9,x) and Pi(9, x) of the reduced phase 
^, defined modulo tt, in terms of the probability distributions Pq{9^x) and Pi{d^x) of 
the accumulated phase 6 e R [15] 

Pii9,x)= J]Pi(^ + n7r,a;) 

Pq{9,x) and P\{9,x) are by definition 7r-periodic, and satisfy the same system (3.7) as 
Po and Pi. As X — > cxD, they converge respectively towards stationary solutions Po{0) 
and Pi (6*) of (3.7) 

Mk - </.osin2^)Po] + noPo(^) - niPi(^) = 

(3.8) 

^ [( A; - (/>o sin 2^)Pi] - noPo (^) + niPi (^) = 
The sum of these two equations gives a simple relation between Po(^) and P\{9) 

{k - (j)osiii20)Po{0) + {k - (j)isin20)Pi{0) = C (3.9) 

where C is a constant which can be evaluated at the point ^ = ^ ; in fact it is exactly 
the density of states N{E) according to (2.18) 

C = kPo + kPi (I) = /cPe, (I) = N{E) (3.10) 

The system (3.8) can now be rewritten as two decoupled equations for Po{0) and Pi{0), 
containing the a priori unknown constant N{E) 

{k-(pi sm20)^[{k - (l)osin20)Po] + {no + ni)[k - Fosm20]Po ^ niN{E) 

(3.11) 

{k-4)osm20)-^[{k-4)i sin26')Pi] + (no + ni)[A; - Pq sin 26']Pi = noN{E) 

Finally N{E) will be determined by the normalization conditions that must be 
imposed on the Ti-periodic solutions Po(6') and Pi{0) of (3.11) 

i:Poie)d0 = ^ ^^^^^ 
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where we have used (3.3). 



According to (2.19) the Lyapunov exponent can be expressed in terms of the 
stationary distributions Po{0) and Pi{0) 

7(£;) =< 0(x)cos2^(x) >= / (i6'cos26'[(/)oPo(6') +</>i^'i(6')] (3.13) 

Jo 



The determination of N{E) and 'y{E) is now reduced to the resolution of the 
equations (3.11) for Po{9) and Po{9). The results are presented in appendix A. In the 
following we will only consider various limiting cases for N{E) and 7(-E). 

* Limit E^0+ 



In order to determine the limiting behaviour of N{E) for £^ ^ 0"*", we have to 
study how the integrals Ii and I2 given in (A. 4) diverge in this limit. It is convenient 



to set fiQ = fxi 



24>o 



and fj, = fii — iJ,Q. 



For fi > 0, Ii has the dominant divergent behaviour and N{E) vanishes as E^ 



N(E) ~ 2a — 
E^o no 



r(^i 



r(/x)r(/.o) 



E 
^0 



(3.14) 



For < 0, I2 has the dominant divergent behaviour and N{E) vanishes as E^^^^ 



N(E) ~ 2(7 — 



r(//o) 



E 
0! 



(3.15) 



• For fj, — Ii and I2 diverge only logarithmically and therefore N{E) vanishes 
very slowly in comparison to the cases 7^ 

N(E) oc (3.16) 

In the next section (IV) it will be shown how these limiting behaviours of NiE) can be 
understood through a qualitative analysis. 
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Let us now turn to the Lyapunov exponent 7(-E'). We can study how the integrals 
Ji and J2 given in (A. 6) diverge in the hmit E ^ 0. Ji (resp. J2) has the dominant 
divergent behaviour if > (resp. < 0) and we recover the result (2.10) 



7(0) = \Fo\ 

The case n — deserves special attention since 7(i?) vanishes as 

-f(E) oc , } (3.17) 

This singularity corresponds to the singularity (3.16) found for the density of states 
N{E) through the Thouless formula [2]. Note that for E = 0, the wavefunctions are 
not exponentially localized since 7(0) = limj^^oo j; In = 0, but it can nevertheless 
be shown that In behaves asymptotically for large L as CL (See Section IV). 

* Limit of white noise process for {0(a;)} 

{(j)(x)} tends to a white noise process | ^ / m ,-,9 c-/ /\ in the 
"-^^ ^ \< (l){x)(j){x') >= F^ + a5{x - x') 

i'o = 00- 



&0 ^ 00 

limit I no 00 provided that 
ni — >• cxD 



V 



^If^ remains constants (3.6). 

a = 2 

no + ni 



The parameter fi = = — also remains constant. In this limit we recover the 

200 cr 
results [8,9] for < < 00 

2a 1 



'^WN{E) = -az£\nCJl{z)+Nl{z) 
where z = ^ and ( J^, A^^) are Bessel functions. 

* Limit i? — > (X) 

In this limit, we recover of course that N{E) behaves asymptotically as the free 
density of states (II 1 b)) 

E 

N(E) ~ 

14 



The Lyapunov exponent ^{E) vanishes asymptotically as ^ in accordance with the 
result (2.16) for exponentially correlated processes. The correspondance between (2.15) 
and (3.5) is the following 



A = no + ni 



and 



noui 



E^oo no + ni \E 



(3.19) 



Let us give a straighforward derivation of this result for the Hamiltonian (See 



Fig.2) 



d 



dx^ 
d^ 



+ 02 + 0'(a;) 

- ^^2+<^O + E(-ir20O<5(^-^n) 

n 

where / ^'^^"^-^ points where 4){x) changes from (— ^o) to (+(/>o) 

1^ {x2n+i) are the points where (/>(x) changes from (+0o) to {—(po) 

On the interval [0, L] there will be typically N attractive delta potentials and N 
repulsive delta potentials with (3.3) 

1 1 

(3.20) 

no ni 

The elementary potentials V± = (j)Q± 20o<^(a^) have respectively the following transmis- 
sion coefficients in the region E > 4>q 



L = N(< l> + <h>) = N(— + — ) 

no ni 



t±iE) 



iCE - (t>l 



(3.21) 



iCE -cPlTcpo 

In the high energy limit, we can neglect interferences between the various delta potentials 
and approximate the total transmission coefficient on [0,L] by a product of elementary 
coefficients t ± (E) [16] 

\t+{E)t-{E)f^^^ (3.22) 



E^oo 



where N{L) 



~ -J Y- according to (3.20). With the help of (3.21) this yields 

L— >oo 1- — 

no ni 



7(E) = — lim Y 111 



^(L) 



^(0) 



1 



E^oo 



_^^\n\t+{E)UE)\ 



(3.23) 



noni 



s-coo no + ni \E 
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Application of the Phase Formalism for E < 



From the spectral properties of H± we know that N{E) identically vanishes for 
E < 0. Nevertheless it is interesting to calculate the Lyapunov exponent ■j{E) in this 
region since it is directly relevant for the problem of classical diffusion [8]. According 
to the dynamical equations (2.21), the Lyapunov exponent ^{E) reads for £^ < 

7(E) = lim ^ / dx[-ksm29(x) + (f){x)cos20{x)] (3.24) 

L^oo L Jq 

and the phase 9{x) evolves according to 

JQ 

— = -/c COS 2e{x) - (l){x) sin 2e{x) (3.25) 



Here, there is no need to distinguish accumulated and reduced phase, according to 
the discussion following (2.21). As before (3.7), we can easily write Master Equations 
for PQ{d,x) and Pi{d,x) 

^ = -M-^ COS 26 - 00 sin 29)Po] - uoPq + n^P^ 

(3.26) 

^ = -^[(-A;cos2^ + 0osin2^)Pi] + noPo - n^P^ 

As a; ^ oo, they converge respectively towards stationary solutions Pq{0) and Pi{0) of 
(3.26) 

^ \{k cos 29 + (t)o sin 20)Po] - n^Po + mPi = 

(3.27) 

^[(A;cos2^-0osin2^)Pi] +noPo -mPi = 
The sum of these two equations gives 

(A;cos26' + 0osin26')Po(6') + (/c cos 26* - 0o sin 2^)Pi(6') = constant (3.28) 

This constant vanishes at the point 6* = f since A;[Po(f ) + Pi(f )] = N{E) = Q. The 
system (3.26) can now easily be decoupled. 

The solutions Po(^) and Pi{9), and the calculation of ^{E) in terms of these 
stationnary distributions are given in Appendix B. Finally, we obtain after some trans- 
formations 

16 



1 / 1 ™ \ '^l 



where <! \ \x + ^y^Jl-x)ix + P^) 



X + 

We have set = ^ , z^i = i^cTTF = WoCTT^- 

• In the white noise hmit (3.6) we recover the result [8] 

K'iz) 



'wn{E = -k') = -2az-p^ (3.30) 



where z = ^ and is the modified Bessel function. 



• In the hmit — > 0, we recover (2.10) 

7(0) = |i^o| 

In section (VI), we will use the expression (3.29) to obtain the velocity for the associated 
problem of classical diffusion in the random force field {(f){x)}. 

IV - LOW^ ENERGY STATES - A PHYSICAL PICTURE 

The purpose of this section is to provide a simple physical picture that accounts 
for the low energy part of the spectrum (3.14, 3.16). We will show that the basic 
mechanism is very different from the one which is responsible for the usual Lifchitz 
singularities. We have to consider separately the cases jji — Q and > 0. 

l)// = 0. 

For a given realization of the disorder, the two independent solutions of iy+V' = 
have been constructed above (2.11). In general, none of these solutions is normalizable 
on the whole line. Instead, if we consider the problem on a finite interval [— i?, i?] with 
Dirichlet boundary conditions, then one can show that there exists a quasi zero mode 
(which is a linear combination of V'o a-nd -01 ) whose energy is exponentially small. One 
obtains with exponential accuracy [19] 



MR) ~ 

it— »00 



1 1 
+ 



rO dx fR dx 
J-RiJ^ Jo 
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^ (4.1) 



Since the corresponding wave function has its support on [—R, i? ] , it is also a 
quasi zero mode on the whole line with same energy. In the presence of disorder, in order 
to estimate Eq wc will replace ipo{x) by its typical value. When {(f){x)} is a random 
telegraph process, the distribution of the random variable dy (j){y) can be obtained 
exactly (see Appendix D). However, for large x, one merely expects that — (j){y)dy 
will be distributed as a Gaussian process [20,21]. More precisely, for the case 0i = — 0o 
and 77-0 = ni = n the Central Limit theorem gives iIj^q^ ~ e~'^°^^~^ . This gives 



x — »oo 



E oc ^e-^-^"^^ (4.2) 

Phrased differently, this means that a quasi zero mode of energy E has typically a 
spatial extension 2R such that 

ln£;~-20oC- (4.3) 

n 

Therefore the number of such states per unit length behaves typically as 

This argument can be easily generalized to account for other types of correlations. For 
instance, if one assumes that the correlation function has a power law behaviour 

<((/>(x)-(/>(y))2>~|x-2/r (a>0) 



(4.5) 



p — y| — > oo 
one can easily show along the same lines that 



N(E) oc — (4.6) 

E^^{-hiE)^ 

This last result is in agreement with the one obtained in [22] in the context of anomalous 
diffusion. 

2) > 0. 

For simplicity we consider the case where 0i = — 0o with ni ^ no. The profile 
of the potential is sketched on Fig. 2. Each time (j){x) jumps from 0o to — 0o and 
then back from — 0o to 0O) there appears a potential well V-{x) = —2(f)oS{x) which is 
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then followed by a potential barrier V^{x) = 2(f)QS{x). Such a dipole like configuration 
can support localized states which account for the low energy behaviour of the density 
of states. Since ni ^ uq the average distance between such configurations is large. 
Therefore one can safely ignore interference effects and consider a single doublet of 
width a. The potential is given by 

V{x) = (j)l- 2(f)oS{x) + 2<j)Q5{x - a) (4.7) 

An elementary calculation shows that this potential can support a bound state provided 



^ - e ^0 (4.8) 



Here we deal with low energy states such that < £^ <^ 0q, this gives 

E ~ ^^e-^^-^o (4.9) 

which is consistent with the previous inequality if a0o ^1- If we would just consider 
the attractive potential well 

V{x) = (j)l- 2cf>oS{x) 

this would lead to a zero energy state ipo{x) oc e~^^^'^° (consistent with supersymmetry). 
It is the coupling of this state with the repulsive potential barrier 2(j)o5{x — a) that 
increases the energy up to £^ = 0Qe~^°''^°. This picture is indeed confirmed by an 
elementary calculation to first order in perturbation theory. Since the width of the 
doublet is distributed according to the distribution 

p{a) = e{a) m e-'^i" (4.10) 

this implies that the number of such states per unit length is 



m — / £j\ 2^o 



N{E)jxe 200 E =1 \ (4.n) 



Since /j, = ^2<i)o° ~ 'Mo ^^^^ result is consistent with eq. (3.14). 



The above discussion clearly illustrates that the basic mechanism which is re- 
sponsible for the low energy behaviour is the existence of low energy states that would 
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have strictly zero energy if one could ignore the couplings between wells. This mecha- 
nism is very different from the one which is at work in the case of Lifchitz singularities. 
In this case, very small energy can occur only if large regions of space are free from 
impurities. This is just the opposite picture of what happens in our case. 

This discussion can be generalized to the case 01 < < 0o with < >^ 0. Assuming 
< 1 01 1 < 00 and no = ni = n, an exact calculation of the density of states gives again 
a power law behaviour 

-{—-—\ 

N(E) a b2 \m mJ . (4.12) 

E — >0 

However the situation is completely different when 0i = 0. In this case, a single 
doublet configuration of length a ^ ^ 

V(x) = (f)l[e{x -a) + 9(-x)] - (l)oS(x) + ^q5{x - a) (4.13) 

supports a low energy state 



E(a) ^ {l-J . (4.14) 



The dependence on the length a is very different from (4.9), and (4.10) gives now a 
Lifchitz singularity 

N{E) oc expf-J^^ (4.15) 



E^O+ V 2^/E 

This result is in agreement with an exact calculation of the density of states that we do 
not reproduce here. 

For < 01 < 005 one can prove that is bounded from below by 0i; N{E) starts 
now with a Lifchitz singularity at 0i 

N{E) = E<^1 

TT m \ (4.16) 



N{E) (X exp - , 



V - COMPUTER SIMULATIONS AND NUMERICAL CALCULATIONS 
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In this section, we consider various probability distributions of the lengths steps 
and discuss the density of states N{E) and inverse localization length 7(-E') as obtained 
from simulations and/or numerical calculations. 



The procedure to obtain the integrated density of states N{E) by simulation for 
the model (III.l) is the following. We generate the random process {0} with the laws 
(3.1) for the lengths of intervals where (f){x) is constant. Then for each energy E, we 
integrate the dynamical equation (2.5) for the total accumulated phase 6 six) 



N{E) is a self averaging quantity and requires therefore only one configuration {0} to 
be computed. We have generated typically 3.10^ intervals, and we have checked with 
a great accuracy that the N{E) obtained was indeed independent of the particular 
realization of the disorder {0}. 

In order to have a better accuracy, we did not integrate the dynamical equation 
for ${x) step by step, but we rather integrated it exactly on each interval where (f){x) is 
constant. For an interval of length £ on which (f){x) = +(f)o, the initial value 9i and the 
final value ^/ of the phase 6{x) are related through 



The integration requires to separate k < (j)Q from k > (See appendix C). 

In Fig. 3, we display N{E) for (pQ = —cpi = 1 and various values of n = tiq = ni. 
Notice that, as expected, N{E) is a continuous function of E. However, it seems that, 
in general, its derivative is not. 



k — (j){x) sva.2dE{x) 



dx 



and we finally compute N{E) with (2.2) 





(5.1) 



For large n (Fig.3a), N{E) is close to the free case ^ except for very small E (3.16). 
This is easily understood if we consider two successive steps of small lengths (~ ^) in 
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(2.5). After two steps the phase increase is 



dOo + ddi ^ {k — (j)Q sin 29)dxo + (A; + 0o sin 26)dx\ k, kdx (5.2) 
This leads to a quasi free density of states. 

In contrast, small values of n (Fig. 3d) lead to 



/>2 

J,2 



- - (5.3) 



i) a plateau N{E) f if E < 
a) a parabolic behaviour N{E) — ^;r^ if E > (j)Q 

The behaviour i) can be explained by noting that, for small n, 6e{L) — 9e{^) ~ iVf ■ 

Therefore after N steps, we get 

NiE) ^ lim — (-n\ = - (5.4) 

^ ^ JV.L^oo 7tL\2 J 2 ^ ^ 

For E > (/)q, one can show (see Appendix C) that 9{L) fa hence we get ii) 

1 - 02 



iV(E) 



T TT 

2 



A quasi- free behaviour develops only for E > (pQ. 

Now we turn to the Lyapunov exponents and consider the ending angles 7r-periodic 
stationary distributions. Wo{0) dd (resp. Wi{d) dd) is the probability for the reduced 
phase to be in the interval [^, 9 + dd] when the system switches from state 0o to state 
01 (resp. from 0i to 0o)- For E > 0, (2.3) and (2.5) lead to 



CE - (po sin 20 



1{E) = --^^^. T'ln 

2(no + ni) y_x V£^- 01 sin 26* 

The stationarity condition implies the relationship 

^ 2 



[Wo{9)-Wi{9)]de (5.5) 



Wo(^)-/ di de'fo{i)Wi{e')5ie-g{i,e')) (5.6) 



(g deduced from (5.1)) and, after some algebra, the differential equations 

[Wo{e){k - 00 sin 2^)] = no{Wi{e) -Wo{e)) 

do 

^[Wi{9){k - 01 sin 2^)] = ni(Wo(^) - W^id)) 
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(5.7) 



Comparison with (3.8) shows that distributions Wo,i(^) and Po,i{d) are the same up to 
a normahzation constant 

Po,i = Wo,i (5.8) 

In particular, we can demonstrate, via (5.8), the equivalence of the two calculations, 
(3.13) and (5.5), of the Lyapunov exponent 7(-E). (5.8) is somewhat unexpected and 
could be possibly related to the Poisson process we use for the step lengths. 

Distributions Wo{9) and Wi{9) are plotted in Figs. 4 and 5 for two values of n. 
Wo{9) (resp. Wi{9)) exhibits a peak at ^ = ckq (resp. at ^ = cui) , the magnitude of 
which is a decreasing function of n. For small n (Fig. 4), Wq is essentially concentrated 
around ckq. On the other hand, for large n (Fig. 5), Wq becomes practically flat. (5.5) 
shows that, in this case, we must get 'j{E) ^ (free case). The small bump at 9 = ai 
(= tto — f when (po = —cpi) is essentially a memory efi'ect due to the preceding step 
{(t){x) — (pi). With the aid of (3.13) and (5.5, 8), we could compute j{E) numerically 
as well as by simulations. Some results are displayed in Fig. 6 (no = ni) and Fig. 7 
(no = 2ni). We clearly see that 7(0) = |Fo| (2.10) and that ^{E), with E fixed, is a 
decreasing function of no (at least for E not too large). This is in agreement with the 
above discussion. The predicted asymptotic behaviour 7(-E') oc E~'^ (3.19) is reached 
for values that largely depend on no and ni. For instance, with no = 1 = 2ni, it is 
practically attained for E>2. But when no = 2ni = 4, values oi E> 20. — 30. are needed. 

Small no,i values deserve a special comment. As already noticed, Wo(W^i) distri- 
bution is peaked at ao{ai). Approximating Po,i(^) by 5-functions, (3.13) readily leads 
to {4>o = -01 ) 

^iE)^Ccl>l-E , E<cl>l 

(5.9) 

1{E)^0 , E><Pl 

This is precisely what we get by means of computer simulations for no = 0.01 = 2ni. 

However, eq.(2.10) 7(0) = |Fo| seems contradictory with (5.9). Investigating the limit 

E —>■ numerically, we could see that very small E values are needed (~ 10""^ — 10~^) 

to detect a decrease of ^{E) towards its limit |Fo|. So, the behaviour of the Lyapunov 

exponent at small energy appears to be far from trivial. 



We now consider step length distributions whose first moment diverges. More 
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precisely, exponential laws fo,i{i) of eq.(3.1) are replaced by distributions who behave 
as 

^-(i+a) £^oo (0 < a < 1) . (5.10) 

a = 1 corresponds to a one-sided Cauchy law. 

Detailed analytical calculations such as those performed in the preceeding sections be- 
come hardly feasible. However, some simple characteristic features appear in the com- 
puter simulations. For instance, when 0o = —4^1, the density of states, N{E), is given 

by 

N{E)^9{E-cf>l)^^^ (5.11) 

TT 

In fact, one can prove that result (5.11) is an exact one. Indeed, since the first moment of 
(5.10) is infinite, the sum L, of N independent positive random variables all distributed 
according to (5.10) behaves, when A?" — > oo, like [23] 

L ~ 7VV« , a<l 

(5.12) 

L NliiN , a = l 
More precisely, jjt/^ (or j^KIy if ck = 1) is distributed according to a Levy stable law. 

Thus, for E < (f)Q, we get 

(5.13) 

Now we examine what happens for E > (jxQ and call -^i, -^2, • • • , -^jv the lengths of the 
successive steps 

L = £i + • ■ • + £jv - iV^/", (a < 1, N ^oo) (5.14) 

According to (C.4), each £i can be written as 

£i = UiT + S£i , 

(5.15) 

0<6£i<T 

and therefore 

L^T{Y^ni) + J2S£i (5.16) 

i i 
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The upper bound ^ S£i < Nt allows to neglect d£i) in L, and we finally get 

i i 

L = (m + • • • + n7v)T + 0(7VV«) 

e{L) = (m + ■ ■ ■ + nAr)7r + 0{N^/") (5.17) 
N{E) = ^ = ^^ E>cf>l. 

Replacing A^^/° by A^lnA^ in (5.14-17), we easily show that the final result (5.17) is 
still valid when a = 1. Moreover, Thouless's formula [2] implies that 

-f{E) = when E xj)^ (5.18) 

Thus, the case a < 1, (po — —<f>i, is like a "free" one, up to the shift of 0q in energy. 

However, simulations with 0o 7^ — 0i reserve some surprise. In Fig. 8, three 
different events drawn with the same law, a = 0.5 are displayed. We now notice that 
N{E) is no longer a self-averaging quantity. In (5.14) splitting L into Lq (and Li), the 
sum of all the steps lengths where (f){x) = (po (or (pi) and observing that Lq and Li are 
both of the same order (N^/°' or A^lnA^), we get, in analogy with (5.14-17) 

^ LoOjE - cPDCE - 0g + Lid{E - cj>l)CE - 0f ^ ^^^^^^ 
7r(Lo + Li) 

in perfect agreement with our simulations. Of course, Lq and Li now strongly depend 
on each drawing. It is interesting to point out that 

i) N{E) - ^ when E ^ 00 

ii) 02 ^ 02 (5 jg) igg^jg ^ e(£-0g^C-E-</.g ^ ^^^^ ^^^.^ self- averaging. 

To conclude this section, we mention that nothing special seems to happen when 
1 < q; < 2. In particular, (5.11, 19) are not observed and the general behaviour of N{E) 
is quite similar to the one we obtained with the exponential law. 

VI - SOME PHYSICAL APPLICATIONS 

In this section we describe some physical problems to which the above model can 
be applied. As mentionned in the introduction, one of the most interesting application 
is the description of classical diffusion in a one dimensional random medium. In the 
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presence of a force (j){x), the probability density P{x t \ y 0) satisfies the Fokker-Planck 
equation 

' 24.(.)P] (6.1) 



dt dx \dx 

supplemented by the initial condition P{x t \ y 0) -^^ S{x — y) 

As we have already explained in the introduction, the transformation 

P{xt\yO) = e-fy '^^"^ '^'Gs{x t | y 0) (6.2) 



(6.3) 



allows to rewrite (6.1) as the Schrodinger problem 

9Gs _ TT r 
Gsixt\y0)^-^^^5ix-y) 

where H+ = — ^ + (j>'^{x) + <j)'{x) is the Hamiltonian whose properties have been 
discussed in the preceding sections. 

^ /"OO 

The Laplace transform Gs{x,y; E) — dt e~^^ Gs{x t \ y 0) satisfies 

Jo 

{E + H+) Gs{x, y; E) = 5{x - y) (6.4) 

/•OO 

The Laplace transform P(x, y;E)= / dt e~^^P{x t \ y 0) therefore reads 

Jo 

P{x, y; E) = Jy '^^'^ < x \ \ y > (6.5) 

ri+ + rj 

When increasing the mean bias < ^ >, it is known that one obtains a succession of 
phases [10] characterized by different anomalous diffusive behaviour. The characteriza- 
tion of these phases requires the knowledge of the velocity and the diffusion constant. 
Here we just consider the velocity given by [8] 

l-= lim <Pix,x;E)> (6.6) 

V E->Q+ 



In order to express V in terms of the localization problem, we consider the average of 
the resolvant at coinciding points 

L— >oo 



hm - / dx<x\ — ^1 x>=-— ^{8} + i7rp{S) 

L— >oo L Jq rlj^ — c — le do 



(6.7) 
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For negative the density of states p{E) vanishes and (6.7) reduces for £■ > to 



1 /"^ , , 1 , d 



The vefocity is therefore given in terms of the left derivative of the Lyapunov exponent 
at zero energy 

1 - d 

— = hm < P(x,x-E) >= hm — ^(-E) (6.9) 

For the model of rectangular barriers of random lengths described in (III.l) we 
have obtained the Lyapunov exponent in the negative energy region in (3.29). 
It follows that the velocity displays a transition at = 1 

V = Q if < 1 

y = MWzl) if >i (6.10) 

In the white noise limit (Eq. 3.6), we recover the result [8] 

y = if \n\<l 

(6.11) 

V = 2a{\n\-1) if |//| > 1 

It is interesting to rederive this result by performing a direct configurational av- 
erage in (6.6). The zero energy Green function at coinciding points P{x,x; 0+) can be 
written for any configuration of the force field {(^(x)} 

P{x, x; 0+) = 2 / d$, e"' L '^(^^'^'^ (6.12) 

J X 

provided one assumes Fq > [8]. Averaging over the disorder {0(x)} gives 

- = 2/ dL<e-'io '^('^)'^"> (6.13) 

y Jo 



For the two-step model (III.l), (f){x) can only take two values (±0o)j therefore 

/ (l){u)du = (f>o A{L) - (f>o[L - A{L)] (6.14) 
Jo 

where A{L) is the random variable that measures the total length in the state 
(f){x) = +00 during the interval [0, L]. 
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The velocity V can thus be expressed in terms of the probabihty density VL{x)dx = 
Pr{x < A{L) < X + dx). One obtains 

1 r°° 

- = 2 / dL< e2^oL-4<^oA(I,) > 

^ Jo 

poo i*L 

= 2 dL e^^°^ / dx e-^'^°''VL{x) . (6.15) 
Jo Jo 

From the expression of Vl{x) given in Appendix D one recover (6.10) after some tedious 
calculations. 

Another interesting application of this model is the fact that it can be used to 
discuss the spectral properties of the two dimensional Euclidean Dirac operator in a 
special type of random magnetic field . Consider the two dimensional Euclidean Dirac 

operator 

il^ = iai (^^^ + iAt^ + ia2 (^-^ + iAx^ (6.16) 
coupled to the external field 

At = fix) 

(6.17) 

Ax = 

If we set ilj{x,t) = expicvt x{x) where x{x) = {v{x),u{x)) is a two component spinor, 
the eigenvalue equation 

= kip . (6.18) 



then reads 



where 



^ — ucf) = kv 

T (6.19) 

J. I. 
vcp = ku 

dx 



(t){x) ^ f{x) +UJ (6.20) 



After decoupling one obtains 

^- ^ ^ x-k\ (6.21) 



H+ 

The knowledge of the density of states of the one dimensional model therefore permits 
to evaluate the density of states in the quasi one dimensional random magnetic field 
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In a previous paper [24] we have shown that the low energy part of p{E) is enhanced 
with respect to the free case. Although the analytical expressions are of course model 
dependent, the enhancement effect is expected to occur on very general grounds. 



The type of disorder that we have introduced in section III can be motivated 
by the following picture. Imagine that the potential At{x) is created by a line of 
dipoles of same strength distributed randomly on the x axis with alternating signs 
+2/0, -2/0, +2/0, -2/0 • • • The resulting potential f{x) = -/o ^e(x - jji) is an en- 

i 

semble of square functions such that each sample function can only take the value ±/o. 
If one assumes further that the probability to find m dipoles on an interval of length x 
is given by the Poisson process 

(nx)'^ 

p(m,x) = 1 — exp(— na;) (6.23) 

ml 

then f(x) is a random telegraph process. After some calculations, very similar to the 
one presented in [24] one obtains 

N{E) ^ - r N{4>o ^ fo + LO ,4>i^-fo+u; div (x ^—^ (6.24) 



(-i"f)' 



Similar logarithmic behaviours have been obtained by a variational method in the con- 
text of the Schwinger model [25]. However a direct comparison with our result is not 
possible because the gauge field measure is not the same. 

Interestingly enough, this two dimensional Dirac operator can also model other 
systems. For instance it arises in the study of the electronic properties of polyacetylene. 
In this context 4'{x) is proportional to the dimerization pattern of the carbon-hydrogen 
chain [26,27]. Introduction of the disorder can of course influence the nature of the 
ground state [28] . It can also model certain layered structures with a piecewise constant 
order parameter. Application of these ideas to Fermi superfluids was recently discussed 
[14]. 

Another potential application of this model is in the field of dynamical systems 
perturbed by noise. A large amount of works have been devoted to the study of two 
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dimensional linear stochastic systems of the form 



x{t) = Ax{t) + r]{t)Bx{t) (6.25) 
where A and B are 2x2 matrices , r]{t) a random noise and x{t) = i ) is a two 



dimensional vector. In order to understand how the linear system x{t) = Ax{t) is 
perturbed by the noise in the asymptotic regime t — > oo one introduces the Lyapunov 
exponent 

A = lim -log||a;(t)|| (6.26) 

t—K30 t 

which characterizes the stability of the solution and the rotation number 

p= hm - arctg (6.27) 

Although there are general theorems regarding the existence of such limits, in most 
cases one cannot compute them analytically. Therefore one has to resort to perturbative 
methods (for instance small noise expansions) [29]. 

The Dirac equation (6.19) can be written in this canonical form with the identifi- 
cation 

^-ce[\ J) B.(; \) (e..B) 

A more physical picture is to consider that x(t) = characterizes the two 

dimensional phase space of an harmonic oscillator perturbated by a random telegraph 
process. This linear system (6.25) can be put into an hamiltonian form with 

H=^{p'' + q'')-mPQ (6-28) 

From the discussion in section III it follows that A and p are nothing but the inverse 
localization length and density of states of the localization problem. 

VII - CONCLUSION 



In this work we have studied a model of localisation based on the Witten Hamilto- 
nian of supersymmetric quantum mechanics. This model exhibits interesting behaviour 
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of the density of states at low energy. We have shown how this behaviour can be un- 
derstood through simple physical arguments. In the case /x = 0, it is the existence of 
quasi-zero modes that is responsible for the logarithmic singularity, whereas for /U ^ 
it is the formation of low energy bound states that explains the power law behaviour. 
These mechanisms rely of course crucially on the supersymmetric structure of the po- 
tential. The resulting expressions are qualitatively similar to the ones that had been 
obtained before in the case of a white noise. In that respect, the existence of a finite 
correlation length does not make much difference. In contrast, when the lengths of the 
barriers are drawn with a distribution with infinite moments, we have demonstrated, 
analytically as well as numerically, that the density of states develops a very different 
behaviour. A similar observation was recently made in the context of anomalous diffu- 
sion. It was shown in [30] that the asymptotic regime can be qualitatively different in 
the presence of correlated transfert rates (sec also [12]). An extension of this model in 
higher dimension is certainly worth pursuing. In particular it would be interesting to 
extend the zero mode analysis to the case of the random magnetic field problem. 



APPENDIX A : CALCULATION OF N(E) AND 7(E) FOR E>0 

The determination of N{E) and ^{E) is reduced to the resolution of the equations 
(3.11) for Po{0) and Po{0), which contain the expressions {k — 0osin2^) and {k — 
sin 29). The solutions will be different if these expressions can vanish or not for some 
values of 9, and it is therefore necessary to separate the cases E > (pQ from E < (pQ and 
E <(j)l from E > 0^. 

In order to simplify the discussion, we will set 0i = —00, in which case the 
spectrum is only divided into two regions < < 0q and E > (f)^. 



We introduce the angle a &]0,^[ such that sin2Q; = CM- and the function 
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cos(6 — a) 




sin(S — a) 


^0 


sin{e+a) 




cos{6-\-a) 





where vi = 



24>Q cos 2a 



and i/Q — 



no 



2(j)o cos 2a ■ 



In this region of the spectrum the resolution of (3.11) requires some care since the 
expressions {k±(pQ sin 2^) can vanish and this is why the yr-periodic solutions Po{0) and 
Pi{0) are defined by two different expressions on the intervals [a, f +«] and [f +«, ir+a] 



a<e<f +a Po(e) = i^N(E) , Jf] , \l 

— —2 ' u\ / 4^,^ \ / cos(a+0) sin(a— J 5- 



^+a<e<.+a Po{e) = fhN{E)- 



{ a<e<^+a P,{e) = f^N{E)- 



f 



^ 2 ' — — ' -1- V / \ / cos(9—a) sin(9+a) Jn — a cos(a+t) sin[a—t)q(t) 



dt 


coa{a—t) 


sin{a+t)q{t) 




dt 


cos(a— t) 


sin(a+t)q'(t) 




dt 


cos(a+t) 


sin{a—t)q{t) 




dt 


cos(a+t) 


sin(a— t)g(t) 



(Al) 



(A2) 



The functions Po{d) and Pi{0) are presented on Figs. 4 and 5 for a certain choice 
of the parameters (f)o, uq, n\. 



The normalization condition J^'^^ Po{9)d0 



N{E) 



N{E) = 



2(7 



— — gives the density of states 

{A.3) 



h + h 

where Ii and I2 are double integrals that read after some transformations 



/i = + / dx 



X 



1 1 r,^fa+/?;r-'|^_,p 



{x + (3'^Y^ Ix-ll'^o {x-l) 



y 



{AA) 



h = (l + r) / dx 



yJ'l 



dy 



{x + P'^yo+^ \x-l\''i Ji y'^i+i y-l 



We have set a = and /3 = tan2Q; = tSP-^- 

no+ni C4>Q—E 



According to (3.13) the Lyapunov exponent ^{E) reads 

/"7r+a 

7(E) = 00 / d^cos2^[Po(^) - i'i(^)] 
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From the expressions (A.l) and (A. 2) for the stationary distributions Po{0) and Pi{0) 
in the interval < < 0q we obtain after some transformations 



^{E) = N{E)[Ji-J2] 



{A.5) 



where 



Ji 



dx 



x^ + /3^{2x-l] 



i: 



x'^ + (3"^ 
1 



1 



- 1 



y 



Vi 1 



x{x-l)y + (3'^ X + (3'^ y{y - 1) 



and J2 is obtained from Ji by the exchange between vi and vq. 



b) Results for E > 4>l 



(^.6) 



1 ^E 

Here we use the parameter — = C—^ > 1 and the function 

9o 



7 



7 ^1 
CI -72 (po 



and 



g(6') = exp ^ — cti arctan( ) — arctan( ^^"^^7 )^ with cti 

CKo = ^ — :tt~ ) ^iid where arctan stands for the principal determination of tan~^ 

CI - Y 00 
taking values in [— f, f ] . 



The solutions Po(^) and Pi(^) of (3.11) read for 9 e [-f , f] 



Po{e)^^N{E)^ 



m 



^ - sin 26' 

7 



(i+sin2t)g(t) 



Pi(^) = ^iV(£;)T 



^ +sin26' 

7 



(i-sin2t)g(t) 



+ B 



where the constants A and B are choosen to ensure the tt periodicity of Po{0) and Pi{0) 



dt 



+ 



1 - e-(«o+«i)'^ (i + sin2t)g(t) e^'^o+ai)^ -1 Jo + sin2t)g(t) 



dt 



and B is obtained from A by replacing + sin2t) by — sin2t) in the integrals. 



TT 

N{E) is obtained by the normalization condition jj^ d9Po{6) — and "fiE) 

by the expression 

7(E) = 00 rf6'cos2^ {Po{e) - Pi{9)) 
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Let us now consider the limit E ^ (f)Q. N{E) is given by two very different 
expressions for E < (j)Q and E > but it is nevertheless continuous at the point 
E = (pQ, as it should by definition of the integrated density of states. N{(f)Q + 0) = 
N{<I>1 - 0) = ^ where 



Ii = Jq dx e ^ / — e y 

J X y 

J TOO , ^-unx r dy -n 

l2= L dx e - ''o^ / e y 
° Jo y^ 



Ml 



{A.8) 



However one should not expect to have more than continuity of N{E) at this point. In 
particular the derivative p{E) = ^ has no reason to be continuous at this point (See 
V). 

APPENDIX B : THE LYAPUNOV EXPONENT FOR E<0 

Here the dynamical equation for the phase 9 is given in (2.21) 

— = -/c COS 2e{x) - (j){x) sin 20{x) (B.l) 
ax 

We have already seen in the discussion following (2.21) that 9{x) is trapped in an interval 
of the form [— |, 0] [modulo tt]. 

For the two step model I '^'^ ^ , I defined in (HI.l), a more accurate discussion shows 
that the phase 6 is trapped in the interval [ai, ckq] [modulo tt] where ckq e] — J, 0[ is 
the angle defined by 

and CKi = — f- — ctQ. 



Indeed on an interval where (pix) = (po the phase 9{x) is attracted towards ao according 
to II = —C(j)Q + sm(29 — 2ao) and on the interval ^(a;) = —(po, the phase 9{x) is 
attracted towards ai according to ^ = —CcpQ + k'^ sin(2^ — 2ai) 



The stationary distributions Po{9) and Pi{9) have therefore their support on the interval 
[ai,ao]. 
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Using the relation (3.28), the system (3.27) can be rewritten as two distinct equa- 
tions for Po{0) and Pi{9), whose solutions are for 6 & [cti, ckq] 



Pl{0) = ^sin(2^e-2ai) 

where 

q{e) = I tan(6' - ao)r° I tan(6' - ai)\''^ 



The constant A is determined by the normalization condition 



(S.3) 



(BA) 



""(i^Po(^)= (S.5) 
ai "-0 + ni 



The expression (3.24) for the Lyapunov exponent ^{E) 

1 

^{E) — lim — / /c sin 26'(x) + 0(a;) cos 26'(x)] 

can now be rewritten in terms of the stationary distributions Po{9) and Pi{9) 

7(£;)= f " de [-ksm2e{Po{e) + Pi{e)) + (f)ocos2e{Po{e)-Pi{e))] 

Jai I J-. 

= C(j)l + /cM de [cos(2^ - 2ao)Po{e) + cos(2^ - 2ai)Pi{e)] 

APPENDIX C : PHASE EVOLUTION ON AN INTERVAL WHERE 0(a;) = 0o 
The integration of (5.1) requires to separate k < (j)Q from k > (j)Q. 
• For k < (po we introduce cto G ]0, ^[ satisfying sin2Q;o = 



dx 



dx 



> <(=^ 6 [modulo tt] e] — ^ — CKo, cko [ 
< <^=^ 9 [modulo tt] e ] ao, f ~ cto [ 



Therefore \df — 9i\ cannot exceed tt. 

More precisely, if {9i + ckq) belongs to the interval ] — f + nir, ^ + mr] with n E N then 
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{$f + CKo) stays within the same interval, and the corresponding tangents are related 
through 

tan(6'/ + do) - tan2cio = [tan(6'j + ao) - tan2Q;o]e"^V'^o-'^ (C.l) 

• For k> (j)Q the increment (^/ — ^j) is not limited. 
The local primitive of the function involved in (5.1) reads 



dO 1^1 tan (/) - 7o \ , 0o 

— arctan — where 7o = ^ 



k-(j)osm20 y/E-(f)l \ ^1 - 7^ / k 

An increment of tt for the phase 6 requires therefore the length r given by 

^ d9 _ 7T 

k-^Q sin 29 ~ CE - 4>l 

If we denote by m G N the integer part of (^), the final phase reads 



^ - I Z J. „:„oz) - 12 (^-2) 



6»j = 6'j + m7r + (/? , (C.3) 
where the angle <^ e [0, tt [ is determined by 

1 



£ = rriT + 



arctan 



tan ^ — 7o 



.-1 



Here arctan is a continuous determination of tan , and therefore 

I tan ^/ — \ 

arctan 



tan — 7o 



APPENDIX D : PROBABILITY DENSITY Vl{x) 



e [ 0, TT [ for e [ 0, TT [ . 



(C.4) 



Let {(j){x)} be a random process which alternates in 2 states (t){x) = (po and 
0(a;) = —00- We denote by 1^ {hn) the successive sojourn lengths in state 0o (resp. 
— 0o)- We are interested in the distribution of the random variable A{L) which measures 
the total length in state on the interval [0, L] (see for instance [31]). The total distance 
in the state (t){x) — —(po on [0, L] is then simply [L - A(L)]. 



The probability density Vl{x) such that 

VL{x)dx = Pr{x < A{L) < x + dx} 
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naturally splits into four terms according to the values of at the end points 0(0) = ±0o 
and (j){L) — ±0o- In obvious notations one has 



VLix) = V++{x)+V+-{x)+VZ+{x)+VZ-{x) (Dl) 

Let us consider the case 0(0) = 0o, (piL) = (po- A typical configuration of this 
type is represented in the following figure. 



The probability density V^~^ (x) therefore reads 
rt^{x)dx = Pr{<j){^) = 0o) X 

oo / N \ ( ^ \ 

^Pr{L-x-dx<^hi<L-x\Pr\'^lj<x<^lA {D2) 

N=0 \ i=l J \j=l j=l J 

The probability involved in this expression can be obtained easily. We have already 
seen that (3.3) 

Prim = = 



no + ni 

Each £i is distributed according to the probability density (3.1) 

Me) = e{l) no e-"°^ 

whose Fourier transform reads 

dx fo{x) e-'P- = 
-oo ^0 + ^P 
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N 



The total length y = is a sum of independent identically distributed random 



i=l 



variables. It is distributed with the probability density 



/o(p) 



N 



= 0{y) no 



N-l 



(N-iy. 



-noy 



{D3) 



(N N+1 \ 

ii < X < ^ ii 1 then reads 
1=1 i=l / 



AT 



Q^ix) = Pr{J2i.<x\-Pr{J2i.<x] = I dy fi^^^y) - fi^+'Hv) 



e{x) 



.1=1 
m 



i=l 



dm. 



N 



,—nox 



We also define fi^^{x) and q\'^\x) by replacing no by ni 

/ N \ 

/i ^''^ (x) is the probability density to have | ^ /i^ = a; j 

(AT AT+l \ 

hi < X < ^ I 
i=l i=l / 



Eq. (D2) becomes 



ni 



no + ni 



Ar=o 



In a similar way we can compute (x) such that 

Vl~{x)dx = Pr{x< A{L) < x + dx \ (^(0) = 0o, = -0o} 



by considering configurations represented in the following figure. 
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We get 



rt~{x)dx = Pr{(f){0) = (f)o) J2 



N=0 

oo 



7V+1 



< ^£i<x + dx 



X _ 

i=l 



Pr 



■ N N+1 

hi < L — X < hi 

.i=l i=l 



Tin + Hi 
^ ^ JV=0 



We can now obtain V]^'^{x) [resp. 7^^ (a;)] from {x) [resp. ■P^"'"(ic)] by the 

simple exchanges ( ^ ^ J^^ so that the probabihty density Vl{x) (Al) finally 
\ X < > L/ — X 

reads 

oo 

'^o + ni L J 



_^ no 



N=0 

oo 



ACKNOWLEDGEMENTS 



We thank Jean-Philippe Bouchaud for many helpful discussions. We also thank Frederic 
Fougere for his collaboration in the early stage of this work, and Olivier Martin for 
reading the manuscript. REFERENCES 

[1] P.W. Anderson, Phys. Rev. 109 (1958) 1492. 

[2] B. Souillard, "Les Houches", XLVI, J. Souletie, J. Vanimenus and R. Stora, 

(1987). 

[3] E. Witten, Nucl. Phys. B188 (1981) 513. 

[4] R. Dutt, A. Khare and U.R Sukhatme , Am. J. Phys. 59 (1991) 723. 

[5] A. Comtet, A.D. Bandrauk and D.K. Campbell , Phys. Lett. 150B (1985) 

159. 



[6] A. Inomata and G. Junker , "Proceedings of International Symposium on 

Advanced Topics in Quantum Physics" (J.Q.Liang, M.L.Wang, S.N.Qiao and 
D.C.Su Eds.) Science Press, Beijing, (1993) 61. 

39 



N.G. Van Kampen, "Stochastic Processes in Physics and Chemistry", North 
Holland, Amsterdam (1981). 

J. P. Bouchaud, A. Comtet, A. Georges and P. Le Doussal , Ann. Phys. 201 
(1990) 285. 

J. P. Bouchaud, A. Comtet, A. Georges and P. Le Doussal , Europhysics. Lett. 
3 (6) (1987) 653. 

B. Derrida and Y. Pomeau , Phys. Rev .Lett. 48 (1982) 627. 

J.P. Bouchaud and A. Georges, Phys. Rep. 195 (1990) 127. 

G. Oshanin, S.F. Burlatsky, M. Moreau and B. Gaveau, Chem. Phys. 178 
(1993). 

C. W. Gardiner, "Handbook of Stochastic Methods", Springer Verlag, New- 
York (1982). 

D. Waxman and K.D. Ivanova-Moser, Ann. Phys. 226 (1993) 271. 

M.M. Benderskii and L.A. Pastur, Sov. Phys. JETP 30 (1970) 158. 

I.M. Lifshits, S.A. Gredeskul and L.A. Pastur, "Introduction to the theory of 
disordered systems", John Wiley and Sons, New York (1987). 

L.A. Pastur and E.P. Feld'man, Sov. Phys. JETP (1974) 241. 

P. Erdos and Z. Domanski, "Analogies in Optics and Micro Electronics" 
(W.van Haeringen and D.Lenstra Eds.) Kluwer Academic Publishers (1990) 
49. 

G. Barton, A.J. Bray and A.J. McKane, Am. J. Phys. 58 (1990) 751. 

40 



E. Tossati, A. Vulpiani and M. Zannetti, Physica A. A164 (1990) 705. 

G. Theodorou and M.H. Cohen, Phys. Rev. B 13 (1976) 4597 ; A. Bovier, 
J.Stat.Phys. 56 (1989) 645. 

J. P. Bouchaud, A. Comtet, A. Georges and P. Le Doussal , J. Physique 48 
(1987) 1445. 

B.V. Gnedenko and A.N. Kolmogorov, "Limit Distributions for Sums of In- 
dependent variables", Addison- Wesley, Reading, MA, 1954. 

A. Comtet, A. Georges and P. Le Doussal, Phys. Lett. B208 (1988) 487. 

A. V. Smilga, Phys. Rev. D46 (1992) 5598. 

H. Takayama,Y.R. Lin-Liu and K. Maki, Phys. Rev. B 21 (1980) 2388. 
D.K. Campbell and A.R. Bishop, Nucl. Phys. B200 (1982) 297. 

B. C. Xu and S.E. Trullinger, Phys. Rev. Lett. 57 (1986) 3113. 
K.A. Loparo and X. Feng, Siam J. Appl. Math 53 (1993) 283. 

C. Aslangul, N. Potier, P. Chvosta and D. Saint-James, Phys. Rev. E47 (1993) 
1610. 

S. Karlin and H.M. Taylor, "A first course in stochastic processes". Academic 
Press, New- York (1975) 



41 



FIGURE CAPTIONS 



Fig. 1 The random force 0(a;) takes alternatively two values 0o and 0i on inter- 
vals whose lengths £'^s and h'^s are independent random variables, distributed 
according to (3.1) unless otherwise stated. 

Fig. 2 The Schrodinger potential (0^ + (j)') as a function of x in the case 0o = — 0i- 
Delta functions appear each time (f){x) jumps. 

Fig. 3 Computer simulations of the average density of states N{E) for 00 = — 0i = 1 
and no = ni = n, with respectively a) n = 5. b) n = 1. c) n = 0.5 d) n = 0.1 

N{E) appears to be an increasing function of n. In general, its derivative is 
not a continuous function of the energy; a discontinuity can occur aX E — (f)^. 

Fig. 4 The 7r-periodic stationary distributions Wo{9) and Wi{9) are plotted, respec- 
tively in full curve and in dashed curve, for the case 0o = ~0i = 1; E = 
0.4, no = ni = 2. Wo and Wi{9) defined just before (5.5) are the same dis- 
tributions as Po{0) and Pi{0) defined just before (3.8) up to a normalization 
constant (5.8). Note the peak of Wq{9) at 6* = ckq ~ 0.34, and the peak of 
Wi{9) at ^ = CKi = CKo — f . For further explanations, see text. 

Fig. 5 The same as Fig. 4 except for no = ni = 4 instead of no = ni = 2. The 
lengths of the intervals on which (f){x) is constant are shorter, and the peaks 
at ao and ai are therefore attenuated. 

Fig. 6 Lyapunov exponents 'j{E) are plotted for (po — —(pi = 1 and uq = rii = n, 
with respectively a) n = 1. b) n = 2. c) n = 4. For these three cases, ^{E) 
vanishes at zero energy, but very slowly (3.17). 

Fig. 7 The same as Fig. 6 except for no = 2ni = n and 7(0) = |Fo| = 1/3, according 
to (2.10). 

Fig. 8 The density of states N{E), for the case 0o = 1 and (pi = — C2, when the 
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steps lengths are drawn according to the same broad distribution, behaving as 
(5.10) with a = 0.5, and whose first moment therfore diverges. Curves a) b) 
c) represent three difi'erent simulation events. Obviously, N[E) is no longer 
self-averaging (5.19). 
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